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1.0  INTRODUCTION 


This  interim  technical  report  describes  work  done  in  year  one  (12  May 
1986  -  11  May  1987)  of  a  program  (Contract  No.  F19628-86- C-0079 )  to 

increase  the  spectral  resolution  of  LOWTRAN  6^^  by  developing  a  radiative 
transfer  model  and  attendant  computer  subroutines  for  molecular  absorption 
effects.  Specifically,  the  goals  for  this  effort  are  to: 

•  develop  algorithms  providing  at  least  3  cm  1 2  resolution  (FWHM 
-  full  width  at  half  maximum), 

•  model  molecular  absorption  of  atmospheric  molecules  as  a 
function  of  temperature  and.  broadening  density  from  0  to 
17,900  cm"1, 

•  treat  the  following  atmospheric  molecules:  water  vapor, 

ozone,  carbon  dioxide,  methane,  nitrous  oxide,  carbon 
monoxide,  molecular  nitrogen  and  oxygen,  and 

•  integrate  the  new  algorithms  and  subroutines  into  LOWTRAN  6. 

The  upper  limit  of  the  spectral  range.  17,900  cm-1,  is  determined  by  the 
HITRAN  line  atlas  for  atmospher  ic  molecules .  ) 


1.  F.  X.  Kneizys,  E.  P.  Shettle,  W.  0.  Gallery,  J.  H.  Chetwynd,  Jr.,  L. 

W.  Abreu,  J.  E.  A.  Selby,  S.  A,  Clough  and  R.  W.  Fenn,  "Atmospheric 
Transmittance/Radiance:  Computer  Code  LOWTRAN  6,"  AFGL-TR-  83-01 87 , 

Air  Force  Geophysics  Laboratory.  Hanscom  AFB,  MA  01731,  (August  1983). 
ADA137786 

2.  R.  A.  McClatchey,  W.  S.  Benedict,  S.  A.  Clough,  I).  E.  Burch,  R.  F. 

Calfee,  K.  Fox,  L.  S.  Rothman,  and  J.  S.  Garing,  "AFCRL  Atmospheric 
Absorption  line  Parameters  Compilations,"  AFCRL -TR  -73-0096 ,  Air  Force 
Geophysics  Laboratory,  Hanscom  AFB,  MA  01731  (January  1973);  L.  S. 
Rothman,  App]  ,  Opt.,  20,  791  (1981).  AD762904 
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A  substantial  fraction  of  this  work  has  boon  completed.  Molecular 
transmittance  through  the  atmosphere  has  been  modeled  using  standard  band 
model  techniques , (35 )  jn  the  present  approach,  the  spectral  region  is 
partitioned  into  1  cm-*  bins  for  each  molecule.  Within  each  bin, 
contributions  from  lines  whose  centers  are  within  the  bin  and  from  nearby 
lines  centered  outside  of  that  bin  are  modeled  separately  with  their  own 
temperature  and  pressure  dependences;  see  Figure  1.  The  absorption  due  to 
lines  within  the  bin  is  calculated  by  numerically  integrating  over  the 
Voigt  line  shaped*.  The  Curtis -Godson ^3 4 5  6  ’ 7 8 9 ^  approximation,  which  is 
accurate  for  the  moderate  temperature,  pressure,  and  concentration 
variations  found  in  the  earth's  atmosphere,  is  used  to  replace  multilayered 
paths  by  an  equivalent  isotropic  one.  The  transmittance  algorithms  have 
been  coded  and  integrated  into  LOWTRAN  6.  Comparisons  to  "exact" 
calculations  with  FASC0D2^*^)  resulted  in  a  reexamination  of  and  a  new 
definition  for  the  standard  line  density  band  model  parameter,  1/d;  this 
new  definition  lowers  the  line  center  absorption  by  taking  the  finite  bin 
widths  into  account.  A  complete  band  model  tape  has  been  generated  for  the 
molecules  listed  above.  A  preliminary  version  of  MODTRAN ,  this  moderate 
resolution  LOWTRAN  6,  has  been  delivered  to  AFGL  for  evaluation. 


3.  D.  Auding,  "Band  Model  Methods  for  Computing  Atmospheric  Slant  Path 
Molecular  Absorption,"  Rpt.  No,  7142-12-T.  Willow  Run  Laboratories, 
Environmental  Research  Institute  of  Michigan,  Ann  Arbor,  MI  (1967). 

4.  A.  J.  LaRocca  and  R.  E.  Turner,  "Atmospheric  Transmittance  and 
Radiance;  Methods  of  Calculation,  Rpt.  No.  107600-10  T,  Environmental 
Research  Institute  of  Michigan.  Ann  Arbor,  MI  (June  1975). 

5.  W.  L.  Wolfe  and  G.  J.  Zissis,  "The  Infrared  Handbook:  Revised 
Edition,"  The  Infrared  Information  and  Analysis  ( I R I A )  Center, 
Environmental  Research  Institute  of  Michigan,  Ann  Arbor,  MI  (1985). 

6.  A.  R.  Curtis,  Q.  J.  R.  Meteor.  Soc ,  78,  165  (1952). 

7.  W.  L.  Godson,  J .  MeteorT,  _1_2.  123  (195*5)  . 

8.  S.  A.  Clough,  "Preliminary  User  Instructions  for  FASC0D2 , ”  Letter, 
AFGL/OPI ,  Hanscom  APB,  MA  01731  (1  October  1985). 

9.  H.  J.  P.  Smith,  D.  J.  Dube,  M.  E.  Gardner,  S.  A.  Clough,  F.  X.  Kneizys 

and  L.  S.  Rothman,  "FASC0DF,  -  Fast  Atmospheric  Signature  Code 
\Spectral  Transmittance  and  Radiance),"  AFGL-TR-78-G081 ,  Air  Force 
Geophysics  Laboratory,  Hanscom  AFB ,  MA  01731  (1978).  APA057506 
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Figure  1.  Spectral  Absorptivity  for  Lorentz  bine  Shapes. 

Curves  Generated  for  Products  of  Line  Strength  and 
Absorber  Amounts  Equal  to  0 . 1  .  1.0  and  10.0  cm~*. 
The  Halfwidth  has  been  fixed  at  0.1  cm"*.  In  the 
Hand  Model  Transmittance  Formulation,  Absorption  Due 
to  Lines  Whose  Centers  Are  Within  a  Bin  (Between  the 
Dashed  Lines)  is  Modeled  Separately  from  That  Due  to 
Line  Tail  Absorption  (Outside  the  Dashed  Lines). 


An  overview  of  the  MODTRAN  code  is  presented  in  Section  2.  In  the 
subsequent  two  sections,  calculation  of  the  band  model  parameters  and  the 
transmittance  formulation  are  described.  Modifications  and  additions  to 
LOWTRAN  6  are  presented  in  Section  5.  In  Section  B ,  MODTRAN  output  is 
analyzed  for  the  purpose  of  validation.  Conclusions  and  future  plans  are 
discussed  in  the  final  seel  Ion. 


The  following  indexing  convention  is  used  throughout  the  report: 

Symbo I  Indexed  Quantity 

i  spectral  bin 

m  molecule 

n  temperature 

cr  species 

I  atmospheric  layer 

The  species  symbol  c,  will  refer  to  either  an  individual  species  such  as 
water  or  ozone,  or  a  mixture  of  molecules  such  as  the  uniformly  mixed  gases 
(UMG)  .  The  atmospheric  layers  4  are  characterized  by  a  fixed  temperature, 
pressure,  and  molecular  constituency.  Constants  used  include: 

c  =  speed  of  light  ( 2 . 997924580xt010  cm  sec  *) 

k  =  Boltzmann  constant  ( 1 . 380662x10 erg  K-1  ) 

h  =  Planck  constant  (6 . 626176xl0-^7  erg-sec) 

N  =  Avogadro  constant  ( 6 . 022045xl0^3  mole'^ ) 

Ts  -  Standard  temperature  used  in  the  HI TRAN  line  atlas  (296  K) 

TQ  =  Standard  temperature  used  in  LOWTRAN  (273.15  K) 
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2.0  THE  MODTRAN  CODE  -  AN  OVERVIEW 


MODTRAN  is  the  name  given  to  the  LOWTRAN  6  code  affixed  with  a  new 
moderate  resolution  radiative  transport  algorithm.  The  new  subroutines 
have  been  written  in  portable  ANSI  I  standard  FORTRAN  and  constructed  so 
that  their  interfacing  with  LOWTRAN  6  minimizes  coding  changes.  These 
additional  elements  do  not  interfere  with  the  regular  operation  of 
LOWTRAN  6;  rather  they  represent  an  additional  capability  for  higher 
spectral  resolution.  The  only  modification  to  the  input  data  sequence  is 


5  -  10  KM  RT  15  0  --  US  Std  Rtm 


2000  2100  2200 


Wavenumber  (cm’*) 


Figure  2.  Transmittance  Curves  for  a  15°  from  Zenith  Slant  Path 
from  5  to  10  km  Through  the  US  Standard  Atmosphere.  The 
Solid  Curve  was  Obtained  from  MODTRAN  by  Passing  a  5 
cm-1  FWHM  Triangular  Slit  Over  the  1  cm-1  Calculation. 
The  Dotted  Curve  is  the  Regular  LOWTRAN  6  Result  also 
Obtained  from  MODTRAN  (20  cm-1  Resolution  with  l)V  Set  to 
5  cm-* ) . 


.in  atlililjnii.il  pnramc  ter .  JHMOl).  which  is  added  to  the  cml  ol  CAKDl .  ’tlic 
user  wliti  is  not  interested  in  ns  i  up.  tin*  MOUTKAN  option  simply  ignores  this 
paramo  1 1-  r  amt  prepares  his  input  file  in  the  usual  way;  the  pi  opi  um  then 
t unctions  in  its  "normal"  mode  and  dot's  not  require  access  to  an  external 
data  fill*  of  hand  model  parameters.  When  usinp  the  hand  model  option 
(.1HMUU  •  1).  artless  to  this  data  file  is  required. 

A  comparison  of  MOUTkAN  and  hOWTHAN  (>  demonstrates  the  importance  of 
the  hip.her  resolution  algorithm.  f'ip.ure  2  shows  transmittance  calculated 
at  h  ami  20  cm  '  resolution  for  a  low  altitude  slant  path  throuj'h  the  US 
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MOUTKAN .  It  the  b  em  '  calculations  were  depended  1 1'  the  20  cm  *  l.OWTKAN  ti 
calculation,  the  overall  upreemeiil  would  be  quite  oil .  However,  at  higher 
resolution,  larpe  deviations  from  the  lower  resolu  »>u  r>l<  illations  occur. 
In  particular,  the  MOUTKAN  calculation  resolves  tb.  structure  of  water 
lines  below  2 HU)  cm  ^  and  tie*  hand  center  of  the  ‘  o  .adamentnl  at  2220 


3.0  M0LK«'.111,AK  HAND  MODEL  PARAMETERS 


The  bus  i  i:  assumption  of  Uu-  bond  mode)  r.iiiu'fi'l  *  ^  is  that  spi'i  I  r.d 

1  i  nr  posit  ions  and  strengths  are  d  i  s  t  r  i  bu  t  ed  in  a  way  that  nan  be 

represent  ed  It  y  a  simple  matliein.it  i  i:a  1  model.  in  mn-  approach.  I  In*  broad 
frocplrnry  c -in',  i  nunni  ul  the  band  models  is  divided  into  small  bins.  '1  be 
line  and  tail  contributions  t  o  each  tun.  Figure  1,  are  then  paramet  i  r  i  /.ed 
tat  th.it  the  bin  t  ransm  i  1 1  ance  ran  he  calculated  when  t  lie  absorber  amounts 
arc  Known 

We  have  chosen  a  bin  width  of  one  wavenumber ,  Av  -  1  cm  ^  .  Line  data 

from  t  lie  111  TUAN  tape^^.  the  Al-'CII.  line  atlas,  is  used  to  calculate  the  hand 

model  paramet  ms  Since  it  has  data  for  lines  in  tin:  t  requeni  y  raitjjc  0  to 
17000  cir  17tP)0  bills  must  he  clmi  art  cr  i/.ed .  Eor  eat  h  emit  i  Unit  i  up, 

molecule,  t  einperat  lire  dependent  absorption  coefficients  and  line  densities 
alonp.  with  a  lint-  width  parameter  are  used  to  paramet  cr  i  /.e  the  spectral 

lints  whose  centers  fall  within  a  ijiven  bin;  a  sini;te  temperature  dependent 
absorption  coefficient  parameter  determines  the  contribution  to  a  bin  f’-om 
the  tails  ol  lines  centered  in  nearby  bins. 

In  the  next  two  subsections,  the  calculation  of  the  molecular  bund 
model  parameters  is  described,  mid  the  f  or  mat  t  i  tit;  of  the  data  file  is 
discussed  in  the  third  subsection.  Formulas  used  to  derive  the  parameters 

are  ijiven.  alonj;  with  a  discussion  of  their  dependence  on  temperature  and 

pressure 

3.1  Line  Center  Parameters 

Each  wavenumber  bin  corresponds  to  a  1  cm  interval  and  contains 

parameters  foi  molecules  with  1  lues  in  tli.it  interval  .  The  molecules  lor 

whirli  hand  tnotlel  parameters  are  included  are. 

11^0.  CO.,.  U, ,  N^t,  CO.  Cll.j.  0^.  tt  N-, . 

The  melerolar  absorption  coeffit  iei.ts  (S/d)  (cm  ^  am. quits  ^  .  1  am.ii'.at 
1  atm  at  SIT')  are  calculated  at  a  reference  temperatures 


IV  ii-1 


v  200  K,  225  K,  250  K.  275  K,  300  K) 


In  the  band  model  subroutines,  a  linear  i nt erpo lot inn  is  used  to  calculate 
absorption  coefficients  at  temperatures  between  200  and  300  K.  For 
temperatures  below  200  and  above  300  K,  the  extreme  values.  (S/d)j  and 
(3/d).  .  respectively,  are  used.  The  molecular  absorption  band  model 
parameters  are  calculated  from  the  individual  line  strengths. 


(S/d), 


Z 


Here  Sj(T,,)  is  the  integrated  line  strenp.th  nt  temperature  T,,  of  the  j’th 
line  of  molecule  m  in  bin  i.  The  line  strength  at  an  arbitrary  temperature 
(Mil  he  scaled  from  the  111THAN  line  strength  at  its  standard  temperature, 


i’ ,  l,y 


Qr(Ts)Qv(Ts) 

Sj(i,,)  =  Qrnn><W 

1  exp(  he V  =/kT  )  V.  ■  T„  T 

sj,t-' 


where  Q,.  and  Qv  are  the  rotational  and  v’hrational  partition  functions,  and 
Kj  is  the  energy  of  the  lower  transition  state. 

A  collision  broadened  or  l.orent/.  line  width  parameter  (tc)y  is  defined 
at  SIT  (T  -  T()  •  273.15  K,  I1  -  1*(|  -  1013.25  mbars).  A  single  val''e  can  be 
stored  because  tin*  pressure  and  temperature  dependence  of  t  he  l.orent/.  line 
width  is  mis i 1 y  mode  1 ed  . 


■yr(  r.n 


‘Vo  7>-  (To/T,? 


whei-e  the  exponent  x(n  has  been  set  to  1/2  loi  all  molecules  except  f'Og,  for 
which  x  •  3/4.  The  (T  .)  band  model  parameter  )s  calculated  as  a  line 
strength  weighted  average  over  the  tabulated  l.orent/,  line  widths  tT,.);(T  ' 


‘Vo  -  <Ts/To>Xm  lX(Vj(Ts)  Sj(Ts)i  /  i^W1  •  (M 

j  j 

i. ike  the  absorption  coefficient  s ,  the  line  density  band  model 
parameters  (l/d)n  (cm)  tire  calculated  at  the  5  reference  temperatures  and 
interpolated  when  used  by  t  tie  band  model  suht  out  i  ties .  The  line  density  is 
del  i tied  by 


(1/d)  = 


N 

V 

>Z  s, 

j  =  l 


7  Z  sj2 


This  definit  ion  fur  the  line  spacing,  which  is  derived  in  the  appendix, 
differs  from  the  usual  definition  involving  a  sum  over  t he  square  root  of 
the  line  si  rengths^'*  •  ^  ^  ^  . 

3 . 2  Line  Tail  Parameters 

The  line  tail  parameters  consist  of  line  contributions  from  lines 
located  outside  of  a  given  bin  but  within  ±  25  cm  '  .  The  1  ‘  ne  tall 
absorption  coefficient  band  model  parameters  C()  (cm  1  umagat  ])  are 
determined  by  Integrating  the  Lorent;;  line  shape  over  ttiis  interval. 


M  Z  (1  Ski’ 

K  -  i  25 


l(S/d)„  Kt-(T„,l»0))k 
(k  i)2  i  1/4 


f I  Ik  i )  Av  | 


where  the  delta  fund  inn  serves  to  exclude  the  k--i  term  from  the  sum  (  i  . »:■ . 
the  line  center  contribution).  For  molecules  other  than  llpl)  and  CCf,  ,  tail 
contributions  beyond  25  cm-'  from  a  line  center  have  been  assumed  to  make  a 
negligible  contribution  to  the  transmittance  and  within  the  25  cm  '  limit  a 


I),  f.  Unhurt  son,  b .  S.  Hernstein,  and  K.  Ilaimes,  "Addition  of  a  5 
cm  1  Spectral  Resolution  baud  Model  Option  to  I.OWTRAN  5,"  AR1  HR  252 
Aerodyne  Research,  Inc.,  hiilerica,  MA  (October  I'JHO). 

1)  0.  Hubert  sou,  I,.  5.  Hernstein,  R.  Ilaimes,  J.  Wunderlich,  and  1, . 

Vega,  "5  cm  ‘  Hand  Model  Option  to  I.OWTRAN  5  ,  App  1  . _ Opd.,  20,  32  Hi 
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form  factor,  f|(k  i)Av],  equal  l.o  1  has  been  assumed.  The  visual  bOWTRAN  6 
water  continuum  consists  of  the  tail  contributions  beyond  25  cm  1  plus 
extrapolated  values  of  the  contribution  within  25  cm'1  (for  smoothness). 
For  CO2 ,  tlie  continuum  from  the  FASC0D21*1  • 11  ^  has  been  added  to  Cn  to 
account  for  the  tail  contributions  from  lines  beyond  25  cm  1 , 


i:n  1  vit:,nh<  2kT-i  — 


Here,  t:  ( v  ■  )  is  the  frequency  interpolated  value  from  FASC0l)2's  block  da  t  a 
/FT.02/.  For  both  llgfl  and  CC^.  the  value  of  C  has  also  been  reduced  by  an 
amount  e<|ua  1  to  the  25  cm  1  from  line  center  contribution  since  this 
contribution  is  included  already  in  the  continuum  data.11'11  The  Cn  are 
taken  to  hr-  proport  ioiuil  to  pressure 


3 . 3  Parameter  Data  File 

Because  of  1  tit'  larp.r  amount  of  data,  the  file  for  the  hand  model 
parameters  is  an  external  file  written  111  binary  format  which  allows  for 
quicker  access  due  in);  the  calculation.  Kacii  entry  corresponds  to  <1  1  cm  1 

interval  and  contains  a  molecular  parameter  set.  Zeroed  data  art*  not 
stored  for  molecules  emit,  r  i  tiu  t  i  nj>  no  lines  to  a  j;  1  veil  interval.  Therefore, 
information  which  labels  the  parameter  set  must  also  be  stored. 

first  entry  of  a  parameter  set  is  the  bin  number,  i.  From  tin*  bin 
number,  the  midpoint  of  the  interval  is  easily  calculated 

Vj  i  Av  (10) 

and  all  lines  whose  centers  fall  in  the  half  opened  interval  ( v  -  Av/2 .  Vj 
1  Av'2j  contribute  to  bin  i. 


12  S.  A.  tllotq;ll.  personal  common  i  cn  V  i  on  . 


A  number  m  between  1  and  8  designates  the  mo  i  cm  1 1* .  Kxcept  for  N2  .  wo 
follow  tile  cmivunl  (oil  set  by  tin?  11  [TRAN  tapn^*^ 

n»  1  2  3  4  b  6  7  H 

molecule  H20  C02  03  N~0  CO  CM4  02  N2 

Tin;  next  entries  in  the  parameter  set  art:  the  molecular  absorption 
coefficients  (S/d)„  (cm  4  amagals  4  )  calculated  at  the  five  reference 
temporal  ui  es .  These  entries  are  followed  by  the  STP  horent/.  half  width, 
(Tc)q,  multiplied  by  1C)4  and  stored  us  an  integer.  bine  spacing  parameters 
(l(j)  for  the  five  reference  temperatures  complete?  the-  line  center  parameter 
set  s . 

bach  parameter  set  for  line  tails  contains  information  on  one  or  two 
molecules.  These  line  tail  parameter  sets  use  the  same  format  as  the  line 
center  parameter  sets.  Attain,  the  first  entry  is  the  bin  number  i  and  the 
second  entry  is  the  molecule  designation  m.  To  recognize,  that  these 
parameter  sets  denote  line  tail  contributions,  their  molecule  labels  are 
offset  by  8 


m 

<) 

10 

11 

12 

13 

Id 

15 

10 

mo  1  ecu l e 

5I20 

C02 

°3 

n2o 

CO 

ch4 

°2 

N2 

coat  i  iniuii. 

parameters , 

<v 

are 

stored 

in 

place 

of  the  (S/d),,. 

Tht 

molecular  designation  and  continuum  parameters  for  a  second  molecule  eau  be 
entered  iu  place  of  the  half  width  ami  (  1  /ct )  parameters. 


4.0  BAND  MODEL  TRANSMITTANCE  FORMULATION 


4.1  L ine  Ce nter  T ransmi ttance 

The  band  model  transmittance  formulation*10*"'  developed  for  the  5 
cm  1  option  to  1.0WTKAN  5*1'1'  has  been  used  to  create  a  model  ate  resolution 
option  lor  1.0WTKAN  (i.*'*  The  expression  used  to  calculate  molecular 
transmittance  due  to  line  centers  is  based  on  a  statistical  mode  i  i«r  a 
1  ini  to  number  of  lines  in  a  spectral  interval,  and  is  given  by 

•t  (1  <Wsl>Mu)<n>  .  (11) 

where  t  is  the  transmittance,  <WS«>  is  the  path  averaged  single  line 
equivalent  width  for  the  line  strength  distribution  in  a  spectral  interval, 
and  <n>  is  the  path  averaged  effect  ive  number  of  lines  in  the  bin 

<n>  --  Au  <  1  /d->  t  i  2 ) 

<  1  / d >  is  the  path  averaged  line  spacing. 

lot  large  <u>  l  (S./d)  and  Av  f  ixcdj  .  the  t  ransm  i  ttaace  simplifies  to  a 
more  reeogn i/.ahlc  exponential  form  Riven  by 

t  ^  exp (  <Wyt>  <l/d>  )  for  <n>  large  (13) 

For  the  molecular  species  iind  relatively  low  temperat  tires  encountered  in 
the  earth's  atmosphere,  the  average  number  of  lines  in  a  bin  is  usually 
small  so  that  the  power  law  transmittance  approximation  is  preferred. 

'lb, ere  are  many  approximations  available  for  calculating  the  equivalent 
width  of  a  Voigt  line  shape;  different  ones  are  valid  fur  different 

13.  F.  X.  Ktieizys,  F.  P.  Shuttle,  W.  0 ,  Gallery  i.  11.  Chetwynd,  .Jr.,  h. 

W.  Ahreu ,  J  F.  A.  Selby,  !<  .  W.  Fenn,  and  i{.  A.  Met'  I  a  t  cliey . 

"  At  mosjihe r  i  e  Transm  i  t  t  am  e/ltad  i  once  :  Gomputer  Code  l.OW'l  <AN  f> .  " 

AFGl,  TU  SO  0007,  Air  Foret;  Gt;opbysics  Laboratory,  Hnnscom  AFB ,  MA 
01731  (February  lOfiO).  ADA08H215 
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regimes.  Doppler  01  collision  broadening,  weak  1  i  in;  or  st  rout',  lin*-,  etc. 

owever,  no  single  approx imat ion  is  adequate  for  t  lit:  range  ot  pressures  and 
optical  path  lengths  encountered  in  atmospheric  transmission  colonial  ions . 
Rather  than  incorporat i ug  different  approximations,  we  directly  evaluate 
the  exact  expression  for  the  equivalent  width  of  a  single  line  with  a  Voigt 
line  shape;  <Wsj>  is  given  by 


<Ws»> 


exp  {  (Su/dJ  -V 1  n2 /tt  K(X.Y}/Ord/d>}dX 


(Mil) 


1-lX.Y) 


Y  f  exp  (  T2)  (IT 
Y2  *  (XT)2 


(14b) 


—  -V  ln2  <n>/<Td/d> 


(140) 


■Win'd  <Tc/d>/<Y(j/d> 


U4d) 


where  F(X,Y)  is  the  Voigt  line  shape  function  and  <Td/d>  and  <Tfc/d>  are 
path  averaged  floppier  and  collision  broadened  line  shape  parameters 
respectively.  The  factor  ol  1/2  in  Kq .  (14c)  corresponds  to  the  half  width 
of  the  spectral  interval- 


4.1.1  Curtis  Godson  Approximation 


averages 


c;i  1  ca  1  ated 


with  the 


Cvtr t  i s  Godson 


approx i mat i on .  ^ ^ ^  This  approx imat ion  replaces  an  inhomogeneous  path 

with  a  homogeneous  one  by  using  average  values  for  the  various  hand  model 


14.  C.  It.  Ludwig,  W.  Malkmus.  J.  1..  Reardon  and  .1 .  A  1,  Thomson, 
"Handbook  of  Infrared  Radiation  from  Combustion  Gases,"  NASA  Report 
SI-  30(H)  (HWO). 
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parameters .  The  Curtis  Godson  approximation  is  very  p  lod  tor  paths  where 
tin-  temperature  or  species  pradienl  is  not  particularly  stoop.  This  is 
certainly  the  case  for  atmospheric  paths  where  the  temperature  variations 
for  arbitrary  paths  fall  within  the  ranee  of  200  to  300  K.  The  total 
optical  depth  is  a  sum  over  contributions  from  the  individual  layers  and  is 
S',  i  veil  by 

ISu/dJ  -  (S/d) j  (An) j  (10) 

i 

where  (Au)j  is  the  incremental  absorber  amount  from  layer  I  and  tS/d),  is 
the  absorption  coefficient  hand  model  parameter  at  the  temperature  dietuted 
by  the  layer-  index,  S.  Note,  frequency  and  species  indices  are  implicitly 
assumed  for  this  and  t he  subsequent  equations. 

The  total  optical  depth  is  used  as  the  weight  inn  function  in 
ralculaiinp  t he  path  averages 

■pd/dj  X  (1/<U|  (s/d)|  lAu)*  (lti) 

i 

J  (s/t,)*  (6u)i  (I7) 

j 


<!/<!> 


<Kr/d>  - 


<Td/d> 


1 

ISu/dj 


(l/d),  (S/d),  (An), 


(18) 


These  hand  model  parameters  were  defined  in  the  previous  section,  and  (Tj)  , 
is  the  usual  Poppler  width  (cm 


<V«  -  (Vi. a  lT> 

v  j  _ _ _ _ _ 

=-  ~  V  2  (  1 p.2  )  N1<T/M 


(10) 


with  M  equal  to  the  molecular  mass  (p/mole)  and  a  the  species  index. 


4.1.2  Uniformly  Mixed  Gases 


Thus  lar,  we  have  implicitly  assumed  that  the  species  o  in  LOWTRAN  are 
simply  molecules  m.  In  actuality.  LOWTRAN  6  assumes  a  uniformly  mixed  i;as 
(UMG)  with  the  molecules  Ct^.  N.,0,  CO,  t'H  j ,  and  02  combined  according  to 
atmospheric  mixing  ratios,  X^  .  ^ '  7>  ^ 

molecule:  C02  N20  CO  Cll^  0-, 

Xm  3.3x10  4  2.«xl0  7  7 . 5x10  ®  1.0x10  (’  2.003x10  1 

Band  model  parameters  for  the  lumuni  formly  mixed  gases,  H ->0  and  0^,  were 
calculated  individually  according  to  the  equations  from  the  subsection  3.1. 
For  the  UMG ,  the  individual  band  model  parameters  were  combined  according 
to  their  ratios 


(S/d)nUMG 


2  X"'.  (S/d)n.m 


(20) 


n/‘DnUMG  -  l  (S/d),/"0  l2  /  ^  l  C  (S/d)n>  m2  /  d/d)n.m  ]  (2.) 

in 


<Vo 


CMC 


1  V” 

liMc;  /  xm.  (^c^o.rn 
m 


(S/d) n 


(22) 


(/  •  J  •—  mm  ■  —  — —  .  .  n  -  - 

nd)nM('  -  ~  1/ 2(  ln2)NkT/M,,M(; 


(23) 


with 


■  UMC 


(S/d)„llMt 


;  2_,  Xm  (^/^^o.m 


( 24  ) 


13.  R.  A.  MrClulcliey,  l< .  W.  Venn,  .1.  1;. .  A,  Selby,  1'.  1  .  Volz,  and  3.  S. 
Gari^g.  "Optical  troperl  i  as  of  the  Atmosphere.”  3rd  ed  .  . 
Al-T.Rl  72  0407,  Air  Force  Cambridge  Research  Laboratories,  Ihmseoni 
AFB ,  MA  10731  (24  August  1972).  A0733075 
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Uric,  thu  subscript  m  indexes  molecules  (skipping  ll20.  O3  and  N2  in  the 
sum),  a  subscript  zero  has  been  added  to  (S/d)  for  the  standard  temperature 
interpolated  absorption  coefficient  ,  and  the  superscript  l)M(i  labels  the 
uniformly  mixed  gases  species.  Again,  we  have  suppressed  usage  of  a  bin 
i udex . 

4.2  Line  Wing  Absorption 

The  power  law  transmittance,  eq .  (11).  takes  into  account  only  lines 

which  originate  within  a  given  spectral  inteival,  and.  for  these  lines, 
only  that  fraction  of  the  line  profile  which  falls  within  the  interval  is 
included  in  the  computation  of  the  equivalent  width.  This  approximation  is 
reasonable  in  the  strongly  absorbing  region  of  a  band;  however,  because  the 
absorptivity  is  expressed  in  terms  of  the  local  1  tne  strength  distribution, 
it  becomes  a  poor  approximation  in  regions  where  the  tail  contributions 
from  nearby  lines  dominate  the  contributions  from  weak  or  nonexistent  lines 
within  a  given  interval  (bin).  This  typically  occurs  in  the  center  and  far 
wings  of  a  band  (i.e..  past  the  band  head) ,  and  also  in  spectral  intervals 
containing  no  lines  which  are  in  the  vicinity  of  strong  lines.  For  these 
situations,  the  local  absorption  is  dominated  by  the  accumulated  tails  of 
the  stronger  lines  originating  outside  the  interval.  The  effect  of  line 
wing  absorption  is  included  in  the  transmittance  by  adjoining  an 
exponent  in 1  term 

t  -  (1  <WS|>/Av)<n>  v  lC,,j  (25) 

where  l  C'u  ]  is  the  total  continuum  optical  depth 

I  <u |  ^  (C)  ,  (An) j  (20) 

* 

The  layer  subscript  5  on  C  labels  both  the  pressure  and  temperature. 


5.0  INTEGRATION  INTO  LOWTRAN  6 


Integration  of  the  MODTRAN  subroutines  into  LOWTRAN  G  has  limn 
accomplished  with  minimal  changes  to  the  original  code.  The  interface 
between  the  regular  l.UWTRAN  <i  ami  the  MGI1TRAN  option  is  made  through  calls 
to  two  subroutines  in  the  LOWTRAN  0  subroutine  TRANS.  A  single  call  to 
subroutine  HMDATA  reads  t  tie  lirst  necessary  wavenumber  block  of  band  model 
parameters  and  calculates  wavenumber  independent,  quantities.  Lor  each 
wavenumber,  calls  to  subroutine  HMOD  are  made  once  fur  initialization  and 
then  additionally  in  die  loop  over  atmospheric  layers  that  calculates  the 
mo  locular  t  ransm i t  t  ance . 

5 . 1  MODTRAN  Subroutines 

The  MDDTRAN  subroutines  are  described  m  this  section.  In  total, 

these  seven  subroutines  contain  less  than  7.r>()  lines  of  code. 

5.1.1  Subrout  i  nt*  HMDATA 

Subroutine  HMDATA  is  called  once  each  calculation  by  subroutine  TRANS. 
This  subroutine  reads  tin*  binary  band  model  tape  header,  makes  the  initial 
call  to  subroutine  ('ONI)KN ,  and  calculates  wavenumber  independent  quantities 
for  subsequent  use  by  subroutine  HMDD. 

After  reading  the  file  header,  HMDATA  checks  the  wavenumber  rang** 

requested  by  the  user  to  see  if  it  falls  within  the  range  in  the  file.  If 
the  requested  interval  is  totally  outside  of  the  band  model  parameter 

range,  tin*  program  stops  with  an  error  message.  However,  if  the  requested 
wavenumber  interval  is  only  partially  outside  of  the  allowed  rang.*;,  the 
program  readjusts  the  upper  and/or  lower  wavenumbers  to  the  file's  values 
and  proceeds  with  the  calculation.  The  ne>:t  read  statement  skips  over 

wavenumber  blocks  in  the  file  until  it  reaches  the  wavenumber  block 

containing  VI,  the  initial  wavenumber. 


If  the  option  1KPT  =  4  is  used,  the  next  calculation  hus  the  same 


atmospheric  and  path  parameters  but  a  new  wavenumber  range.  In  this  case, 
the  call  to  BMDATA  only  rewinds  the  file  and  performs  the  initial  reads. 
The  common  block  /CAKD5/  was  added  to  pass  the  value  of  1RPT  to  BMDATA. 

When  using  the  MODTRaN  option,  the  array  WPAT1I  contains  the 
Incremental  absorber  amounts  for  each  layer.  However,  WPATH  is  not 
calculated  by  subroutine  PATH  when  the  path  lies  completely  within  one 
atmospheric  layer,  so  it  mast  be  specified  in  subroutine  BMDATA  for  use  by 
subroutine  BMOD .  Finally,  several  wavenumber  independent  quantities  for 
each  layer  are  calculated  in  BMDATA  and  stored  in  the  appropriate  arrays 
for  later  use  by  BMOU .  These  are  Vt  .  the  pressure  normalized  to  one 
atmosphere  (P/P0),  and  the  temperature  indices  for  interpolation  of  the 
band  model  parameters. 

5  1.2  Subroutine  CONUKN 

Subroutine  CONDEN  makes  binary  file  reads  of  the  band  model  parameter 
tape.  The  tape  contains  information  on  each  molecule  separately.  To  form 
the  uMC  parameters ,  as  required  by  the  L.DWTRAN  6  structure,  subroutine 
CONDEN  also  condenses  the  molecular  information 

Data  is  only  stored  for  molecules  contributing  to  a  spectral  interval; 
therefore,  it  is  convenient  to  use  recurrence  formulas  to  sum  over 
molecules  for  defining  the  HMG  band  model  parameters: 

(S/d)^™*1)  =  <S/d)n<m>  ♦  <s/d)n  m  +  1  (27) 


<l/d)n‘»*l> 


(S/d)n<",1>  ]2  /  {  1  ( s/d)n<m>  ]2  /  (i/d)n<m) 


(28) 


-  IS  - 


( ),,(n"  1 ) 


(S/d),  ,("• >> 


{rrl;)„(m)  (s/d)u 


(in) 


'  Xm '  1  ^c'o.ni'l  (S/cl)o,mi  1  : 


(2<l) 


M  (  m  M  ) 


(  S/.i )  ^  < m  '  1  ) 


(M 


(in ) 


(S/d)„<m)  -  Xnll  ,  Mlni,  (S/d)„  nill) 


(SO) 


wIhtc  the  supersi  r  i  i>t  s  denote  Du-  level  (it  rn-uiTi,m”,1  i.e.  (in)  and  (ni'1) 
aic  tin1  old  and  new  value,  respect  i  vc  1  y ,  of  tin?  (IMS  hand  model  paramet  ers . 
Note,  these  retentions  are  derived  from  Kqs .  (20)  through  (2.4). 

3.1.3  Suhroul  i  ne  PMOl) 


Subroutine  HM01)  calculates  the  transmittance  using  the  statistical 
hand  model  described  in  Section  4.  Ait  initialization  call  is  made  cure 
every  wavenumber  from  subroutine  TUANS,  and  subsequent  calls  ere  made  for 
each  layer  to  calculate  the  transmittance.  On  the  first  ca’i,  r.'.uii  calls 
CONllhN  il  the  next  blocK  of  baud  model  poraiia.t c rs  needs  to  be  read  I1W  is 
the  count  jug  variable)  and  then  zeros  out  quantities  for  via?  transmittance 
calculations.  On  the  subsequent  ca 1 1 ( s )  ,  IlMhl)  calculates  the  total  optical 
depth  and  transmittance  for  each  species.  The  continuum  cunt <  ilmtiuns  due 
to  tails  of  lines  original  inn  outside  the  1  cm  1  interval  are  included  in 
the  calculation.  I1M01)  returns  the  transmittance  for  the  molecular  species 
to  TRANS.  The  t  ransm  i  1 1  anees  are  stored  in  tin*  regular  I.OWTRAN  ti 

transmittance  array  (TX). 

The  call(s)  made  by  TUANS  to  HMOD  for  ca  1  cu  1  at  i lig  the  transmittance  is 
embedded  in  a  loop  over  atmospheric  layers.  When  only  t  rausin  i  t  t  ,iu<  e  is 
calculated,  the  ho  loop  in  TRANS  consists  of  one  pass.  When  atmospheric 

radiation  is  calculated,  the  TK.»NS  loop  is  over  all  layers,  because  the 
radiation  depends  an  tlu:  increment  at  change  in  transmittance  for  each 
layer.  The  Curtis  Godson  approximation  for  the  statistical  hand  model  also 
requires  the  coat  r  Unit  ion  fmm  each  layer  for  transmittance  calculations. 
In  keeping  with  the  philosophy  of  minimizing  changes  to  the  basic  l.OW'i'KAN  li 

program.  this  layer  loop  is  done  in  HMOD  for  transmittance  only 


is  done  in  HMUh  fur 


i;il  culat  inns  and  in  TRANS  lor  radiation  cu  Icul  at  ions .  Thus,  the  layer  loop 
in  HMOD  consists  of  just  unc  pass  when  "1KANS  is  looping,  over  all  layers  ami 
vice  versa.  Subroutine  UMOP  determines  whether  the  loop  in  TRANS  is  over 
oik:  layer  or  all  layers  raid  then  adjusts  its  own  loop  variable  accordingly. 

The  absorber  densities  (Wi'ATH)  calculated  in  1.0WTKAN  are  converted  to 
the  units  of  i'ib  amapats  in  BMOD  to  be  consistent  with  the  units  of  the  hand 
model  parameters.  The  absorber  density  conversion  factors  from  the  repular 
bOWTKAN  b  units  to  those  used  in  MOPTRAN  are: 

SrtTlKS  l.OWTRAN  t>  MODTRAN 


H-,o 

1  p./cm2 

- 

i.ar>  x  to 3  cm 

amapat 

HMD 

1  km 

1.0  x  1  O'*  cm 

amapa  t 

lb, 

1  cm  amapnt 

- 

1  cm  amapat 

wlii'iT  1  amar.at  1  atm  at  STP . 

The  curve  of  growth  used  in  the  statistical  model  is  based  on  the 
equivalent  width  of  a  single  average  line  in  an  interval.  The  equivalent 
width  and  line  center  t  ransm  i  1 1  mires  are  calculated  by  suhrnui  i  ne  i'.W  i  ii'i'ii 
us  i  up,  a  Vuiet  line  shape  to  combine  the  Doppler  and  Loreutz  line  shapes, 
'fhe  Curtis  (aidsiiii  approximation  is  used  to  calculate  the  l.orentz  halt  width 
(SI),  Doppler  half  width  (S2)  and  line  density  (SB)  for  the  equivalent 
homopeneuus  path.  Ttie  effective  nuinher  ot  lines  is  also  calculated. 

In  addition  to  the  molecular  contributions  from  lines  within  each  1 
rni  1  spectral  interval,  contributions  from  tails  of  lines  external  to  that 
interval  are  a  1  so  i  uli  iihdi'il  Since  these  tails  have  a  smooth,  spectral 
structure,  they  form  e.  emit  inuuni  component,  amt  the  simple  exponent  i  a  1  form 
is  used  to  calculate  the  transmittance.  These  tail  contributions  and  the 
molecular  components  are  combined  (multiplied  together)  and  stored  in  the 
mitrix  TN  for  subsequent  use  by  subroutine  TRANS. 


f> .  1  .  4  Subroutine  CALC 


CALC  is  culled  by  subroutine  BMOU  to  determine  t  lie  Imlul  model 

parameters  tor  arbitrary  temperatures.  The  tabulated  band  model  parameters 
are  linearly  interpolated  over  temporal  lire  tor  each  layer.  If  the 

atmospheric  temperature  is  outside  of  the  temperature  ranee  of  the 

tabulated  parameters,  the  entry  for  either  the  lowest  or  highest 
temperature  is  returned. 

b.l.fi  Subroutine  KWIPTH 

KWII'TH  is  called  by  subroutine  tlMOlt  to  ealoulaie  the  equivalent  width 
of  a  single  average  line  for  the  1  em  1  interval  in  addition  to  the  line 
center  transmittance.  As  discussed  in  Subsection  3.1,  the  Voigt  line  shape 
is  numerically  integrated  over  the  full  interval.  When  the  optical  depth 
at  tin*  line  center  is  less  than  XMIN./ir,  Peer's  law  is  used.  Otherwise 
1, WIDTH  compares  the  relative  values  of  the  Doppler  and  l.orenty.  half  widths 
and  selects  the  proper  formula  for  the  numorirai  iiuep.i  at  ion. 

Several  numerical  approximations  are  used  in  cal cul at i ug  the 

equivalent  width.  The  Voigt  line  shape  function  is  the  real  part  of  the 
complex  error  function^*’)  and  is  calculated  in  subroutine  t’i’1'12.  The 
region  out  to  three  half  widths  is  calculated  by  using  a  linear 
approx  i  me,  t  i  on  in  the  integration  interval.  From  a  study  of  various 
combinations  for  the  number  of  integration  intervals  for  each  half  width 
tlSTF.!’)  and  tin*  manlier  of  halt  wi'  t  hs  from  the  line  center  (NAl.F), 
integrating  out  to  three  half  widths  in  six  steps  was  found  to  represent  a 
reasonable  tradeoff  between  numerical  accuracy  and  execution  time.  The 
region  beyond  three  half  widths  of  the  line  center  is  calculated  using  an 
usymptot ic  expression. 


IP.  M.  Abramowity.  and  1.  A.  Stegun,  "Handbook  of  Mai  hemal  leal  Fund  ions 
with  Formulas,  Graphs.  and  Mathematical  Tables."  NPS  Applied 
Matbeniat  irs  Series  ST>  (December  5  )  . 
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5  .  1  .  ft  Subrout  i  no  KRl-'U 


l-.RR!  is  used  by  KtolbTH  to  calculate  tin;  error  (’unction.  A  polynomial 
approximation  which  is  accurate  to  better  than  2.5  x  10  5  is  used.^lfi^ 

5.1.7  Sub' out  i  ne  CIM'12 

01*1*12  is  called  by  EWIDTW  to  calculate  the  complex  error  function^  ^ 
and  its  first  derivative.  The  real  part  is  the  Voij*l  line  shape  function, 
and  the  first  derivative  is  used  to  approximate  its  frequency  dependence 
when  per  form  i  up,  the  integration  in  MWIDTII.  The  complex  error  function 
subroutine,  which  was  developed  by  J.  Humlicek,^^  is  based  on  a  12'th 
order  Hermit e  polynomial  approximation.  An  3'th  order  approximation  was 
investigated,  hut  the  results  were  not  sufficiently  accurate. 


5  2  Mor.ii  i cat  ions  to  LOWTRAN  6 

An  attempt  has  been  made  to  minimize  modifications  to  LOWTRAN  6.  As 
tneut  ioned  in  the  hep,  i  an  i  up,  of  this  section,  the  switch  .1HM0D  has  been  added 
to  /CAHD1/.  Only  if  JRMOI)  equals  1  are  any  changes  initiated. 

A  few  routines  have  required  minor  modifications.  In  the  MAIN,  JHM(H) 
is  read  and  the  bard  model  parameter  file  is  opened.  Also.  MAIN  has  been 
modi!  ied  so  that  a  divide  by  zero  error  does  not  result  when  AhAMl  is 
determined  at  a  frequency  of  0  cm  The  molecular  column  density  in 

5  i  uroi  u,  a  no  rut*  pressure  i  in  e  a  t  u  a i mu  spue  r  i  r .  ]  a  y  e  f  in  G  F.O  a  i‘  i*  now  s  t  ri  rod  in 
common  blocks  lur  use  by  the  MODTKAN  subroutines. 

The  st  ra  ip.lit  forward  application  of  the  hand  model  transmittance 
tormulat  ion  to  the  single  scattering,  model  in  LOWTRAN  t>  would  requ  i  re 
doubling  of  fill  layer-  dimensions  and  storage  of  too  miiny  path  riat;i  points*  . 


IV.  K  •  M.  Goody,  'Atmospheric  Radiation,"  Uxloed  University  Press  (IPM). 
Id.  .1.  limn  1  i  eek ,  "An  Efficient  Method  for  Kva  lur:  t  ion  of  t  he  romp  l  ex 
Probability  function:  The  Voipt  Function  and  i»s  Derivatives." 

.1  .  Uu. nit  .  Spec  t  .  Rad  . _ Transfer ,  21,  909  (197H). 
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This  problem  was  avoided  by  defining  a  single  layer  for  the  primary  solar 
path.  The  single  layer  approximation,  which  is  only  used  by  the  MODTKAN 
subroutines,  requires  the  calculation  of  Curtis  (nutson  averaged  pressures 
and  temperatures  in  SSGKO .  Significant  modifications  were  required  for 
TKANS .  200  lines  of  code  have  been  added  to  this  subroutine,  but  its  basic 
logic  remains  unaltered.  When  JUMOD  is  set  equal  to  one,  TKANS 

•  sets  frequency  step  size  to  1  cm  1  . 

•  bypasses  ClIiTA,  C2DTA.  and  U311TA  and  calls  the  MODTKAN 
subrout l nes , 

•  interpolates  t ransm i l tances  calculated  in  0  cm ' ^  steps  for 
species  other  than  l^O.  CMC.,  and  Oi,  and 

•  employs  a  discretized  triangular  slit  function  with  FWHM  set 
to  INT(l)V)  to  automatically  degrade  the  1  cm'  1  hin  results. 

In  Subsections  5.1.1  and  5.1.3,  the  calls  to  MODTKAN  subroutines  were 
discussed.  Note  that  the  discretized  triangular  slit  function  when  DV  is 
set  to  1  cm  1  is  a  1  cm  1  rectangular  slit  function,  as  shown  in  figure  3. 


FWHM  CONTINUOUS  DISCRETIZED 

TRIANGULAR  TRIANGULAR 

SLIT  SLIT 


Figure  3.  Comparison  of  Continuous  and  Discretized  Triangula*- 
Slit  Functions.  FASC0D2  Uses  the  Weighting  from  the 
Cold  unions  Function,  blit  MODTKAN  Uses  the  Discretized 
Approx i ma  t  ion. 


6.0  VALIDATION 


Validation  of  the  liand  model  imraniot  ers  on  the  data  tape  has  been 
completed.  However,  only  a  limited  amount  of  testing,  of  MODTKAN's 
t  ransm  i  t  t  ane.e  formulation  has  been  pel  formed  thus  far.  All  validation  ot' 
the  code  will  be  by  comparison  to  KASC0D2  results  since  it  is  the  infinite 
resnlut  ion  lint?  by  line  calculations  of  KAK('.ni)2  wbieb  MODTRAN  is  at  t  empt  i  up, 
t  o  approx i ma t  e . 


0-1  Band  Mode:  1  Parameters 


The  new  band  model  parametois  were  validated  apainst  the  moderate 
resolution  bOWTRAN  5  (10, 11)  values.  The  1  cm  ^  square  slit  hand  model 
parameters  ran  be  depended  to  .1  cm  '  square  slit  band  model  parameters  (J 
odd)  using  the  equations 


j  I  (S/.H, 


(31) 
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where  J  Mins  from  it  (.1  l)/2  to  i.li(.)  l)/2.  Tile  labels  i  and  j  lu-re  denote 
hin  i  ml  lies  and  tin-  swpersc  r  i  p  t  s  (1)  and  (.1)  refer  to  the  1  rm  '  and  .1  1'  in  ^ 
ha  lid  mo  ile  1  parameters,  respectively  (Note  that  t  lie  line  tail  absorption 
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Figure  4.  Comparison  of  Band  Mode)  LOWTRAN  5  (Dotted  Curve) 
and  Band  Model  LOWTRAN  6  (Solid  Curve)  Molecular 
Absorption  Coefficients  lor  Water  Vapor  at  300  K. 

coefficient  formula  is  approximate;  au  exact  expression  for  its  spectral 
degradation  is  not  derivable.)  In  Figure  4,  the  spectral  structure  of  the 
absorption  coefficients  is  compare:!.  Using  these  formulas,  the  1  cm  1  band 
mode)  parameters  were  degraded  to  5  ciiT  *  and  then  compared  to  the  band 
model  parameters  used  in  the  moderate  resolution  option  to  LOWTRAN  S . 
Degradation  of  the  new  parameters  to  T>  cm  ~  *  yields  complete  agreement. 

6 , 2  Comparison  with  FASC0D2 

In  Figure  5.  MODTRAN  and  FA3C0D2  transmittance  curves  are  compared  for 
a  low  altitude  slant  path  through  the  US  Standard  Atmosphere  In  the 
spectral  region  2H>t>  2  200  on ' ' ,  The  curves  were  generated  with  a  1  cm  ^ 
rectangular  slit  (non  overlapping  for  MODTRAN  and  overlapping  for  FASC0D2). 


H2u  Vapor  at  300K 


2/BB  ?/?0  2/40  ?/M)  V/m  2tS0B 

WnVENUMBER  (cm*1) 


5  -  10  KM  RT  15° 


US  St d  Rtro 


rsJ 


Wavenumber  (crn  U 

Figure  r> .  Oonipar  i  son  of  a  FASC0D2  and  u  MODTRAN  Calculation.  N^O 
aid  CO2  are  the  Primary  Absorbers. 

The  absorption  in  this  region  is  primarily  due  to  a  N2O  band  centered 
around  2224  cm  '  and  the  P  branch  of  the  4.3  pm  CO-,  band.  Although  the 
agreement  is  very  good,  discrepancies  do  exist.  The  dip  at  2210  em  ' 
stands  out  .  The  N-,0  lines  in  t  li  i  s  region  are  spaced  about  0.8  to  0.0  cm  * 
apart,  and  therefore  most  1  cm  *  bios  contain  just  one  line.  The  Inn  a1 
2211)  cni  ’  is  an  except  inti  cunt  a  i  11 1  ng  two  lines,  resulting  in  the  dip. 
Similar  analysis  can  ho  applied  to  other  discrepancies  and  implies  that 
better  agreement  results  with  some  minimal  degradation  of  the  spectral 
rosu 1  at i on 

Since  the  hand  model  parameters  were  generated  down  to  0  cm  *  ,  MODTRAN 
can  he  used  to  perform  molecu'ar  absorption  ea  lculat  ions  in  the  microwave. 
Results  for  a  0.1  Km  se.i  level  path  through  the  US  Standard  Atmosphere  are 
compared  to  FASCUD2  in  Figure  (i .  The  comparison  is  good,  but  problems  arc 


T ransm i t 1  once 


evident.  A  series  of  lines  at  7 ,  11,  13,  ami  15  em  are  missing  in  tin 

MODTRAN  on  1  tail  at  i  on .  Further  investigation  is  necessary  to  dot  ermi  no  the 


origin  of  these  discrepancies 


0.1  KM  PATH  ot  0  Kf1  ALT  —  US  Std  Rtm 


<w 


Wavenumber  (cm’1) 


Figure  6.  Comparison  of  a  FASC0D2  and  a  MODTRAN  Calculation  in 
the  Microwave. 


7.0  CONCLUSIONS  AND  FUTURE  WORK 


A  radiative  transfer  model  applicable  to  upper  altitudes  lias  been 
developed  1  Hr  1.0WTRAN  C> .  This  moderate  spectral  resolution  model  has  been 
coded  and  integrated  into  LOWTRAN  6.  Testing  has  begun  on  the  new  code, 
referred  to  as  MODTUAN,  and  a  preliminary  version  has  been  delivered  to  the 
Air  Force  Geophysics  Laboratory. 

In  the  following  year,  most  of  our  effort  will  concentrate  on 
validation  and  opt  imiz.it : on  of  the  code.  For  the  frequency  ranee  from  0  to 
17,5)00  cnT 1  ,  MODTUAN  and  FASCOD'g  calculations  will  be  compared  for  low  and 
high  altitude  paths.  Significant  discrepancies  will  be  identified  and 
corrected.  After  validation  is  completed,  the  MODTUAN  code  will  be  made 
more  efficient.  Methods  which  may  be  fruitful  for  optimizing  the  code 
i  nr.  1  ude 

•  defining,  cut  off  values  for  absorption 

coef  f ie i out  s , 

•  combining.  hand  model  paramot  ers  of  absorbing 

species  before  determining  transmittance  rathei 
than  determining  t  rnnsni :  1 1  mice  of  each  species 
separat c ly , 

•  developing  a  more  efficient  routine  for  integrating 
the  Voigt  line  shape,  and/or 

•  allowing  for  the  spectral  resolution  of  the  baud 

model  parameters  to  be  degraded  internally. 

To  the  extent  time  permits,  each  ot  these  approaches  will  be  investigated. 
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APPENDIX:  DERIVATION  OF  THE  BAND  MODEL  PARAMETERS 


In  MODTUAN ,  the  molecular  trnnsmi  tt  mice ,  t,  from  lines  originating 
within  a  spectral  bin,  Av,  is  determined  from  an  expression  of  the  form 


Av/2 

(  dv  )" 

0 


(A  1) 


where  b'v)  is  a  lint-  shape  function,  u  is  the  absorber  amount  ,  and  S  and  n 
are  functions  of  the  absorption  coefficient  (S/d)  and  line  density  (1/(1) 
hand  model  parameters 


(S/d) 

(A  2) 

(1/d) 

«•  (1/d)  Av. 

(A-  3) 

Eq.  (A  1)  results  from  Kqs .  (11)  and  (14a). 

A  standard  method  for  determining  bund  model  parameters  was  laid  out 
by  Goody  ^  ^ .  The  weak  and  strong  line  limits  of  a  two  parameter 
transmittance  expression  are  determined.  These  results  are  equated  to  the 
weak  and  strong  line  limits  of  the  I.ndeuburg  Reiche  function,  and  the 
resulting  equations  are  solved  for  the  unknown  parameters.  When  Goody 
applied  his  approach  to  two  trial  transmittance  formulas,  both  gave 
expressions  tor  the  band  model  parameters  of  the  form 


(S/d) 


N 


(A  4) 


A  1 .  K .  M .  Goody 
( 1964) . 
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(A  5) 
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/ 1  Si 
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Here,  N  is  the  number  of  lines  and  S-  is  the  integrated  line  strength  of 
line  i  (the  line  spacing  formula  assumes  all  N  lines  have  the  same  half 
width) . 

It  is  important  to  note  that  Goody's  approach  is  inappropriate  for  the 
transmittance  expression  in  MODTRAN .  Eq .  ( A - 1 )  models  the  transmittance 
through  a  set  of  lines  whose  centers  fall  within  a  given  finite  spectral 
bin.  The  badeitburg  -  Re  i  che  function,  on  the  other  hand,  models  the  total 
absorption  due  to  those  lines  -  not  just  their  contribution  within  a  finite 
spectral  bin. 

Instead  of  equating  the  weak  and  strong  line  limit  of  Eq.  (A  1)  to  the 
hadenburg  Reiche  function,  we  equate  the  weak  and  strong  line  limit  of  Eq. 
(A  1)  to  (Mass'  expression^  ^)  for  the  transmission  due  to  lines  whose 
centers  are  randomly  distributed  in  a  given  spectral  interval, 

Av/2  Av/2 

2  r  2  r 

(  t,xbl  Snb(v)]  dv  )"  it  Ilj  {  —  I  exp  l  -S  jtibl  v )  1  riv  }  (A  fi) 
0  0 


Here,  all  N  lines  are  assumed  to  have  the  same  half  wiith,  but  the  line 
shape  function  is  arbitrary. 

It  !i  •  i  !•.!»!«  !»''  t  lin  wnt 


Av/2 

a  r* 

~J  l  b(V)]J/j!  dv. 
0 


(A  7) 


then,  in  the  weak  line  limit,  Kij.  (A  6)  becomes 


A  2.  ti .  N.  (Mass,  "Models  for  Spectral  Hand  Absorption,"  ,1.  Opt  Soc . 
Am. ,  48,  GOO  ( (908) . 


1  n  (I  (Rt  II  S)  I  u2  l  ~Y~  (  »  s  )Z  «  (R2  -~2~)  "  ^  I  '  01  )  ~ 

N  B  2  N  2  N 

1  i  u  <  u2  l  (2si)2  +  iu2  -1T)  2  Sj2  1  '  0(l,3)-  tA  n) 

i  1  i -1  » -1 

Tin;  first  order  equation  in  u  reestablishes  Kq .  (A  4) 


N 

"  S  =  I  Si 

i  1 


(A  ',)) 


To  dot  ermine  the  strong  lino  limit  o  f  Kq.  (A  G).  wo  assume  b(v)  to  bo 
a  non  i  noroas  inj;  bounded  positive  function  on  1 0 .  a> )  .  Consider  a  general 
integral  of  the  form 


Av/2 

2L-  C  e-Sltb(V)  Jy 

Al'J 

0 


(A  10) 


It  follows  that 

0  <  ,,  -Sub  (0 )  <  j  <  Snl.(6v/2)  <  j  (A  11) 

This  inequality  requires  that  1  falls  off  exponentially  with  large  u  and 
that  the  exponent  is  linear  in  u.  The  exponential  decay  factor,  0,  is  given 
by  the  limit  as  u  -*  ®  of  the  negative  logarithm  of  1  divided  by  u 

C  ■=  1  im  llnll)/uj  (A  1  * 

Applying  l,' liopilal  1  s  rule,  one  obtains 


C  *  S  <h(v)> 


(A  Id) 


win;  re 


Av/2  Av/2 

b  ( v )  e"Sub(v)  dv  /  J  e“Sub(V)  dv.  (A  14) 

0  0 

Note  that  this  lin-.it  is  independent  of  the  value  of  S  as  long  as  S  >  0. 

Taking  the  logarithm  of  both  sides  cf  Eq .  (A  6).  dividing  by  u,  and 
uquut i m*  tin;  asymptotic  limits,  one  obtains 

N 

v* 

n  S  <b| u ) >  “  ^  si  <b(v)>.  (A  15) 

i  -1 

Since  <li(v)>  is  line  independent,  the  surprising  result  is  that  Eq .  (A  9), 
the  weak  line  limit,  is  also  the  strong  line  limit! 

To  determine  n  and  S.  we  return  to  the  weak  line  limit,  Eq.  (A  ft),  and 
consider  the  equation  obtained  from  the  coefficient  of  the  quadratic  terms 
in  u.  This  equation  gives 
N 

n  S2  =  2^  Si2-  (A  16) 

i  =  1 


<b(  v)>  =  lim 


U-K» 


« 
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Tie;  resulting  expressions  for  tht;  band  model  parameters  are; 

N 

(S/d)  =  T7.  y  si  (A  4) 

nr 

t  - 1 


(1/d) 
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i  -l 
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N 

V 

z 

i  =1 


(A  17) 


(S/d)  has  the  siime  form  as  tie!  ore.  but  a  smaller  value  fur  (1/d)  is 
obtained,  because  the  weaker  lines  are  weighted  less  in  Eq .  (A  17)  than  in 
Eq.  (A  r>)  . 
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Figures  5  and  A-l  eompure  a  1  cm' 3  non  overlapping  rectangular  slit 
MODTRAN  calculation  with  the  new  and  old  band  model  parameters  to  a  1  cm  1 
overlapping  rectangular  slit  FASC0D2  calculation.  The  results  with  the  new 
formulation  are  much  improved.  Further  investigation  and  validation  are 
planned . 
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Figure  A-l.  Comparison  of  FASC0D2  and  MODTRAN  Calculations. 

The  MODTRAN  Curve  was  Generated  witli  the  Old  Band 
Model  Parameters .  Eds  ( A  -  'l )  A  (A-5),  and  no  N2O 
Continuum  Function 
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